Formation of deeply bound ultracold Sr2 molecules by photoassociation near the 

^S + '^Pi intercombination line 
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We predict feasibility of the photoassociative formation of Sr2 molecules in arbitrary vibrational 
levels of the electronic ground state based on state-of-the-art ab initio calculations. Key is the 
strong spin-orbit interaction between the c IIu, A^E^ and B^Si" states. It creates not only an 
effective dipole moment allowing free-to-bound transitions near the ^S -I- ^Pi intercombination line 
but also facilitates bound-to-bound transitions via resonantly coupled excited state levels to deeply 
bound levels of the ground X^E^ potential, with v" as low as v" — 6. The spin-orbit interaction 
is responsible for both optical pathways. Therefore, those excited state levels that have the largest 
bound-to-bound transition moments to deeply bound ground state levels also exhibit a sufficient 
photoassociation probability, comparable to that of the lowest weakly bound excited state level 
previously observed by Zelevinsky et al. [Phys. Rev. Lett. 96, 203201 (2006)]. Our study paves 
the way for an efficient photoassociative production of Sr2 molecules in ground state levels suitable 
for experiments testing the electron-to-proton mass ratio. 
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I. INTRODUCTION 



The cooling and trapping of alkaline-earth metals and 
systems with similar electronic structure have attracted 
significant attention over the last decade. The interest 
in ultracold gases of alkaline-earth atoms was triggered 
by the quest for new optical frequency standards [l|. 
The extremely narrow linewidth of the intercombination 
"'^S -I- ^Pi transition, together with the magic wavelength 
of an optical lattice [2,] , is at the heart of the clock pro- 
posals. Strontium is the atomic species of choice in many 
current clock experiments 3-6]. The narrow width of the 
intercombination line implies Doppler temperatures as 
low as 0.5 /iK for laser cooling [7|]. It also allows for easy 
optical control of the atom-atom interactions via optical 
Feshbach resonances that involve only small losses @, Q . 

The diatomic strontium molecule represents a candi- 
date for high-precision spectroscopy that aims at deter- 
minin g th e time- variation of the electron-to-proton mass 
ratio [101 . The idea is to prepare tightly confined Sr2 
molecules in their electronic ground state by photoasso- 
ciation in an optical lattice and carry out high-precision 
Raman spectroscopy on the ground state vibrational level 
spacings [10, lli| • Photoassociation refers to the excita- 
tion of colliding atom pairs into bound levels of an elec- 
tronically excited state [I4I . Molecules in their electronic 
ground state are obtained by spontaneous decay [l3| . 
Whether the excited state molecules redissociate or decay 
into bound ground state levels is determined by the shape 
of the excited state potential curve and possibly its cou- 
pling to other excited states. Long-range potential wells 
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and strong spin-orbit interaction in the excited state of 
alkali dimers were found to yield sig nificant bound-to- 
bound transition matrix elements [ij]. 

To date, Sr2 molecules in their excited state have been 
formed by photoassociation, using both a dipole-allowed 
transition |15l . Il6j and a dipole-forbidden transition near 
the ^S-l-^Pi intercombination line [13, El- The for- 
mation of Sr2 molecules in their electronic ground state 
has not yet been demonstrated except for the very last 
bound level [Taj ■ After photoassociation using the dipole- 
allowed transition, the majority of the excited state 
molecules redissociates, and only the last two bound lev- 
els of the electronic ground state can be populated [131 ■ 
This is due to the long-range R~^ nature of the electroni- 
cally excited state (with R denoting the interatomic sepa- 
ration) that does not provide any mechanism for efficient 
stabilization to bound ground state levels [20] . The situ- 
ation changes for photoassociation near the intercombi- 
nation line where the excited state potential curve in the 
asymptotic region behaves predominantly as R"^ with a 
small SC^^R'^ correction, where 5C^^^ is proportional 
to a^ (with a the fine structure constant) . Large bound- 
to-bound transition matrix elements with the electronic 
ground state that behaves asymptotically as R~^ are then 
expected [13] ■ However, quantitative estimates on which 
ground state levels can be accessed were hampered to 
date due to lack of reliable ab initio information on the 
excited state potential energy curves and, importantly, 
the spin-orbit interaction. The latter is crucial because 
it yields the effective dipole moment that is utilized in 
the photoassociation transition and also governs possible 
bound-to-bound transitions following the photoassocia- 
tion. 

Here, we consider the photoassociation process of two 
ultracold strontium atoms into the manifold of the cou- 
pled c3n„(iS + 3p) + AiS+(iS + iD) + BiS+(iS + ip) 



states. The excited state potential energy curves, spin- 
orbit coupling and transition dipole matrix elements are 
obtained by state-of-the-art ab initio calculations [2ll |. 
This allows us to make quantitative predictions on the 
photoassociation rates, bound-to-bound transition ma- 
trix elements, and spontaneous emission coefficients. We 
find that the spin-orbit interaction alters parts of the 
excited state vibrational spectrum qualitatively, opening 
the way for transitions into deeply bound ground state 
levels. This implies that the standard picture of pure 
Franck-Condon type transitions near the classical turning 
points in the ground and a single excited state potential 
energy curve yields qualitatively wrong predictions. The 
crossing between the c3nu(iS + ^P) and AiS+(iS + ^D) 
states is found to also significantly affect the transition 
moments for the Raman spectroscopy envisioned for the 
test of the electron-to-proton mass ratio. The paper is or- 
ganized as follows: Section llll introduces our model and 
briefly reviews the theoretical methods employed. The 
numerical results are presented in Section Hill and Sec- 
tion |IV] concludes our paper. 



II. THEORY 

When a pair of colliding atoms absorbs a photon, it 
undergoes a transition from the scattering continuum of 
the X^E+ ground electronic state into a bound rovibra- 
tional level of an electronically excited state. Here, we 
consider photoassociation using a continuous-wave laser 
that is red-detuned with respect to the '^Pi intercombina- 
tion line of strontium. This transition is dipole-forbidden 
in the nonrelativistic approximation. The c'^II state, cor- 
relating to the asymptote of the intercombination line 
transition, is, however, coupled by the spin-orbit interac- 
tion to two singlet states, A-^E+ and B^E+. Both singlet 
states are connected by a dipole-allowed transition to the 
ground electronic state, X^E+. Thus an effective tran- 
sition matrix element is created which for moderate and 
large interatomic separations is well approximated by 



dso 



_ (Xil]+|d,|BiS+)(Bil]+|Hso|c3n„) 

(X^S]+|d,|AiS]+)(Ail]+|Hso|c3n„) 
A . 



E, 



'c^n^ 



- Ef.1 



'AiSi 



(1) 



where Hso is the spin-orbit Hamiltonian in the Breit- 
Pauli approximation j22l ]. The long-range part of dso, 
dominated by the first term in the above expression, is 
due to the coupling with the B^I]+ state, ideally suited 
for photoassociation. The short-range part is due to the 



coupling with the A^S+ state, paving the way toward 
efficient stabilization of the photoassociated molecules to 
the electronic ground state, as we will show below. The 
scheme for photoassociation into the lowest manifold of 
Hund's case (c) 0+ states is depicted in Fig. [TJ 

We will make use of non-adiabatic effects caused by 
the spin-orbit interaction and therefore employ the dia- 
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FIG. 1: (Color online) Proposed scheme for the production 
of ultracold Sr2 molecules by photoassociation near the inter- 
combination line. The green wavefunction represents a scat- 
tering state of two Sr atoms and the red, blue and brown 
wavefunctions the diabatic components of the excited state 
vibrational level with binding energy E^i^^i^ — 12.9 cm~^. 
Spin-orbit interaction facilitates a transition from this level 
to X^S+ v" = 6 (with the corresponding wavefunction de- 
picted in purple) via spontaneous or stimulated emission. 
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FIG. 2: (Color online) Spin-orbit couplings (left) and transi- 
tion dipole moments (right) between the relevant electronic 
states of the Sr2 dimer that enter the Hamiltonian ((2|). 



batic (Hund's case (a)) picture for our calculations. The 
corresponding Hamiltonian in the rotating-wave approx- 
imation reads 
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where dz{n <— X) denotes the z component of the elec- 
tronic transition dipole moment from the X electronic 
ground state to an electronically excited state n. R is 
the interatomic separation. The peak amplitude and the 



J 



detuning of the photoassociation laser with respect to the 
intcrcombination line are represented by Eq and Aj^^ , re- 
spectively. The diagonal terms for the (n)'^^'^+^^|A| state 
are given by: 
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J(J + 1) + 5(5 + 1) - n^ - S^ + L{L + 1) - A^ 
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with ^ denoting the reduced mass, Vn 



|A| 



(i?) the ra- 



dial potential energy curve, J the rotational quantum 
number, and S the electronic spin quantum number. A, 
E, and Q. denote the projections of the electronic or- 
bital angular momentum, electronic spin angular momen- 
tum, and the total angular momentum on the molecular 
axis, respectively. The term involving the electronic or- 
bital quantum number L in Eq. ([3]) is an approxima- 
tion to the true diagonal adiabatic correction [231, with 
L corresponding to the orbital quantum number in the 
separated-atom limit, cf. the discussion following Eq. 
(40) in Ref. [2j]. The spin-orbit matrix elements are de- 
fined by 



A{R) 



(4) 



(c3n„(I] = ±1, A = Tl)|Hso|c'n„(S = ±1, A = Tl)), 

and 

a(i?) = (c3n„(E = ±l,A = Tl)|Hso|AiE+), (5) 
UR) = (c3n„(E = ±l,A = Tl)|Hso|BiE+). (6) 

The potential energy curve for the X^E+ ground elec- 
tronic state was taken from Ref. [2J|. All other poten- 
tial energy curves, spin-orbit coupling matrix elements 
(shown in the left panel of Fig. [5]) , and electronic tran- 
sition dipole moments, dz{n ^ X) (shown in the right 
panel of Fig. [2|), were obtained from state-of-the-art ah 
initio electronic structure calculations. The details of 
these calculations as well as their agreement with the 
most recent experimental data [23, in particular for the 
crucial A-'^E+ state, are reported elsewhere [21j . 



The most promising route to form Sr2 molecules in 
their electronic ground state via photoassociation and 
subsequent spontaneous emission is determined by diag- 
onalization of the full Hamiltonian ([2]) and analysis of its 
rovibrational structure. In order to connect our model to 
experimental observables, we calculate the photoassocia- 
tion rate, Kiuji, T), and the branching ratios for sponta- 
neous emission, P{v" ■<— u'). The absorption coefficient 
K{uji,T) at laser frequency Wi is given by |26l.l27| 

V'J' J" 
/•OO 

X / e-^/^--^|5„v'(£^,^",wi)Pd£;, (7) 

where p denotes the gas number density, T the temper- 
ature, ks the Boltzmann constant, v' and J' the vibra- 
tional and rotational quantum numbers in the electroni- 
cally excited state, J" the rotational quantum number of 
the initial scattering state, gjn the spin statistical weight 
depending on the nuclear spin, equal to one for ^*Sr, and 
Q^ = {p,kBT/2TTh^f/'^. S^'j'{E,J",u}i) is the 5-matrix 
element for the transition from a continuum state with 
scattering energy E and rotational quantum number J" 
into the bound level \v' , J'). Throughout this paper, the 
quantum numbers J" and v" denote the rovibrational 
levels of the ground electronic state, while J' , v' refer 
to the rovibrational levels of the excited electronic state. 
The square of the S matrix element in Eq. ([7]) can be 
approximated by the resonant scattering expression for 
an isolated resonance V2m. 



\S,,j.{E,J",uJi)\ 



r.'.AE,J"H>r 



[E - /\,.j,{uj,)Y + WyAe, J") + t,j,] 



(8) 



where 'y^,j,{E,J") is the stimulated emission rate, 
j'/;,{E,J") the rate of the spontaneous decay, both in 
units of h, Ay>ji{u!i) is the detuning relative to the 
position of the bound rovibrational level \v',J'), i.e., 
Ayiji = Eyiji — fkoi, where Eyiji is the binding energy 
of the level \v\J'). In Eq. (|5]), we assume the decay rate 
due to any other undetected processes to be negligible. 
The spontaneous emission rates 7^/j/ are obtained 



J 



from the Einstein coefficients Ay' 



J',v"J" 



v"J" 



Ay'J\y"J" 



(9) 



and related to the natural lifetimes Tyij', jyiji = h/Tyiji 
The Einstein coefficient Ay'j\y"j" is given by 



A.y> 



J',V"J" 



4a3 



Hj,{Ey'j, - Ey''jMY^{x^„j„\d.{n' ^ ^)\x':'r) 



(10) 



r 



where Hji is the so-called Honl-London factor equal to 
(J' + l)/(2 J' + 1) for J' = J" -1 and J'/{2J' + 1) for 
J' — J" + 1, and e denotes the electron charge. The 
label n' represents all considered (singlet) dissociation 
limits of the excited diatomic molecule, in our case these 
are ^S + ^P and ^S + ^D. The non-adiabatic rovibra- 
tional wave functions x"j(i?) = {R\Xvj) ^^^ obtained as 
the eigenfunctions of the coupled-channel Hamiltonian, 
Eq. (I2|), in the absence of the photoassociation laser field, 
i.e., for Eq — 0. In principle, in Hund's case (a), the rovi- 
brational wave functions x"j(-R) could also be labeled, in 
addition to n, v and J, by the quantum numbers p, 5, 



S, A and 51, denoting the parity, total electronic spin, its 
projection on the molecular axis, the projection of the 
orbital electronic angular momentum and projection of 
the total electronic angular momentum on the molecu- 
lar axis [231 ■ Since here we consider bosonic ^*Sr atoms 
which are photoassociated to form molecules in the rovi- 
brational states of the 0+ potential, the parity is equal 
to one, and the projection of the total electronic angular 
f2' is zero, which in turn implies A'=0 for singlet excited 
states n'. 

At low laser intensity, /, the stimulated emission rate 
is given by Fermi's golden rule expression: 



Yy,j,iE,r)^4n''- E E K* 



M" 



-J" M' 
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\d-€]^y,rM')\^ 



(11) 



where e denotes the vector of the laser polarization, c is 
the speed of light, '$ej"M" and \E'«'j'M' denote the total 
non-adiabatic (electronic and rovibrational) wave func- 
tions of the initial and final states, respectively. M is the 
quantum number of the projection of the total angular 

momentum J on the space-fixed Z axis, and d denotes 
the electric dipole moment operator in the space-fixed 
coordinate system. After introducing the Born-Huang 
expansion of the non-adiabatic wave functions, Eq. pTJ) 



can further be simplified to the following form [2E 



i:,.AE,r) =aa-{2j' + i)Hj 



Y^{x\jAd.Xn' ^X)\xij,) 



,(12) 



where Xej" (R) a-re energy normalized continuum wave 
functions of the ground electronic state with scattering 
energy E. Using this notation, the transition matrix ele- 
ments between coupled-channel rovibrational eigenstates 
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FIG. 3: (Color online) Population of the c^n„, A^E„ and 
B^Eu components of the Oj rovibrational levels for J' — 1. 
The binding energies are taken with respect to the Sr( Pi) + 
Sr(^S) asymptote. 
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{v'\r\d.\v\J') ^J2(^l,j„\d.in' ^X)\x::j,). (13) 
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FIG. 4: (Color online) Vibrational wavefunctions of the cou- 
pled c"^II„, A^E„ and B^Ej electronic states for v' — —6, 
w' = — 15 and v' — —26. The corresponding binding energies 
are S^/^.g = 0.27 cm~^ = 8.09 GHz, E^,^_is = 12.9 cm"\ 
and _B„/=_26 ~ 75.8 cm^^. Note the different scale for the 
interatomic separation. 



They are almost J— independent as a result of the ex- 
tremely small spacings between the rotational levels 
of Sr2. We may therefore assmxie {v" ,J"\dz\v' ,J') w 
{v"\dz\v') (of course, the selection rule J" = J'±l holds). 



Finally, the branching ratio. 



Piv" 



v'J') = 
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(14) 



describes the probability for the spontaneous decay from 
the level \v' ,J') of the electronically excited state to rovi- 
brational levels |w", J" = J' ± 1) of the ground electronic 
state. Again, the branching ratio P{v" <— v' J') is nearly 
independent of the J' quantum number. 



III. NUMERICAL RESULTS AND DISCUSSION 

We consider ®*Sr atoms trapped at a temperature of 
T ~ 2 ^K, typical for the two-color mangeto-optical traps 
employed for the alkaline-earth species [29|. At such a 
low temperature, the collisions are purely s-wave, i.e., 
J" = 0. The Hamiltonian ^ is represented on a Fourier 
grid with an adaptive step size |30l432l |. 

The photoassociation yield is determined by the 
ground state scattering length and the rovibrational 



structure of the levels in the excited c^Il,,, A 






B^E 



state manifold which couple to J symmetry. The correct 
ground state scattering properties including the scatter- 
ing length are accounted for by employing the empiri- 
cal X^S+ potential reported in Ref. [2J|, reflecting the 
current spectroscopic accuracy. The excited state rovi- 
brational levels are obtained from diagonalization of the 
Hamiltonian ([2]) with Eq = 0. Their analysis reveals a 



significant singlet-triplet mixing, cf. Fig. |3] presenting the 
c'^n,j, A^E+ and B^S+ diabatic components of the cou- 
pled wavefunctions. This mixing results from the cross- 
ing between the c'^n„ and A^I]+ states, which are cou- 
pled by spin-orbit interaction. On average, the rovibra- 
tional levels are predominantly of triplet character as ex- 
pected for the ^S H- '^P asymptote. However, a sequence 
of peaks indicates occurrence of rovibrational levels with 
very strong singlet-triplet mixing. These levels are par- 
ticularly useful for both photoassociation and a subse- 
quent bound-to-bound transition. This is illustrated in 
Fig. m showing the vibrational wavefunctions that cor- 
respond to the two right-most peaks in the A-'^S+ state 
components of Fig. [3] (at binding energies of 12.9 cm"-'^ 
and 75.8 cm~^) and comparing them to the v' — —6 wave- 
function, the lowest level previously observed experimen- 
tally [17]. The v' = —6 wavefunction is almost purely 
long-range and of predominantly triplet character, with 
the population of both singlet components being three or- 
ders of magnitude smaller than the triplet one (note that 
the wavefunctions of both the A^S]+ and B^I]+ compo- 
nents were scaled up by a factor of 100 to be visible in 
the figure). The picture changes completely for the lev- 
els v' = —15 and v' ~ —26. Since the relative weights 
of the c^n„ and A^S]+ components are almost equal, cf. 
Fig.[3l the v' — —15 and v' — —26 wavefunctions in Fig. [4] 
display A^I]+ and c'^H^, components on the same scale. 
Remarkably, the triplet wavefunctions also show peaks at 
short internuclear distance. This is a clear signature of 
resonant, non-adiabatic coupling between vibrational lev- 
els of the spin-orbit coupled electronic states [Ij, |3j, ISJ] . 
It occurs when two potential energy curves that are cou- 
pled cross and the energies of the two corresponding 



vibrational ladders coincide [33|. Then the vibrational 
wavefunctions reflect the turning points of the two po- 
tentials, as seen in Fig. HI Resonant coupling was shown 
to lead to significantly enlarged bound-to-bound transi- 
tion rates to form deeply bound molecules in their elec- 
tronic ground state [ij, IH, [S^l . According to Fig. SI it is 
the coupling between the c^Hu state and tH"e"'A-'^E+ state 
that becomes resonant, inducing strong mixing between 
these components. The effect of this resonant coupling 
will be further increased by the presence of the B^E+ 
state in addition to the A-'^S+ state. The behaviour of 
the B^E+ component strictly follows the c'^n„ wavefunc- 
tion, but is two orders of magnitude smaller, cf. Fig. ID 
This is easily rationalized in terms of the B^E+ com- 
ponent representing only a small admixture, due to the 
spin-orbit coupling ^2(^) in the Hamiltonian 1^, to the 
principal part of the (1)0+ state that originates from the 
c^n„ potential. The magnitude of the B^I]+ component 
is straightforwardly estimated by treating the spin-orbit 
coupling as a perturbation and calculating the first-order 
correction to the wave function, similarly to the expres- 
sion for the transition dipole moment, Eq. ([1]). 

In the alkali dimers, the spin-orbit coupling mixes in 
a triplet component that does not directly participate in 
the optical transition between singlet states 14, 35, 36]. 
Therefore, the enhancement of the bound-to-bound tran- 
sitions in the alkali dimers is only due to the modi- 
fication of the singlet wavefunction. Here, for bound- 
to-bound transitions to the electronic ground state, the 
effective dipole is mainly due to the coupling between 
the c^n„ and the A^J:+ states, cf. Eq. ([T]). There- 
fore, it is not only the modification of the c^Hu wave- 
function but also the presence of a large A^E+ compo- 
nent that is responsible for the enhancement of bound-to- 
bound transitions. Both effects together, the additional 
peaks in the c^n„ wavefunction at interatomic separa- 
tions R < lObohr, and the large A^E+ component at 
these interatomic separations lead to a significantly en- 
hanced effective dipole moment according to Eq. ^ . We 
thus find that for alkaline-earth atoms near the ^S -I- ^Pi 
intcrcombination line, the resonant coupling enlarges the 
singlet admixture to a predominantly triplet wavefunc- 
tion and enhances both the bound-to-bound and the free- 
to-bound transition matrix elements. The enhancement 
of the bound-to-bound transitions significantly reduces 
the lifetime of the excited state bound levels. The life- 
times of the levels v' = —15 and v' — —26 are found 
to be 30.9 ns and 27.2 ns, respectively, compared to 7.61 
/iS for v' = —6, i.e., they are decreased by two orders 
of magnitude. This is rationalized by a larger sponta- 
neous emission rate resulting from an enhancement in 
the bound-to-bound transitions according to Eq. ([TU)) . 

The two effects, i.e., an increase in the bound-to- 
bound and free-to-bound transition matrix elements, 
have an opposite impact on the photoassociation prob- 
ability, with the former hindering and the latter facil- 
itating the photoassociation process. The photoassoci- 
ation absorption coefficient, cf. Eq. ([7]), is shown in 
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FIG. 5: (Color online) Photoassociation into rovibrational 
levels of the coupled c^n„, A^Ej and B^Ej states below 
the Sr(^Pi) + Sr(^S) dissociation limit for a laser intensity 
1 = 1 W/cm^ and two temperatures, 2 fiK (upper panel) and 
20 ^K (lower panel) . The transitions to the six least-bound 
levels that were reported in Ref. [17t | are shown in a semi- 
logarithmic plot in the insert (note the different scales). 



Fig. [S] for all bound levels below the ^S -I- ^Pi dissoci- 
ation limit for two temperatures, T = 2 yuK [29| and 
T = 20 fiK [141 . The absorption coefficient for the levels 
that were experimentally observed [17] are shown in the 
inset of Fig. [5] using a logarithmic scale. At T = 2 /iK, 
the peak rate coefficients for the strongly mixed levels 
v' = -15 and v' = -26 amount to K = 1.6 x lO'^^ 
cm'^s^^ and K = 1.8 x 10^^^ cm'^s""'^, respectively, com- 
pared to K — 1.9 X 10"^^ cm'^s^^ for the lowest pre- 
viously observed level, v' — —6, i.e., about one order 
of magnitude smaller. However, at T = 20 //K and also 
at higher temperatures, the levels with strong resonant 
coupling have absorption coefficients that are very sim- 
ilar to that of v' = —6, K — 1.3 X lO^^^cm-^s^^ and 
K = 1.6 X lO-i^cm^s-i for v' = -15 and v' = -26, 
respectively, compared to K — 2.2 x 10""'^^ cm^s~^ for 
v' — —6, see also bottom panel of Fig. [5l The peak 
rate coefficients for the strongly mixed levels are less af- 
fected by temperature broadening. This is rationalized 
in terms of their large natural width, of the order of a 
few MHz. In constrast, for the level v' = —6 the natu- 
ral width amounts to merely 20 kHz. The natural widths 
need to be compared to thermal widths of 42 kHz and 
0.42MHz for T = 2/zK and T = 20 /iK, respectively. 
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FIG. 6: (Color online) Shape of the photoassociation line for 
transition into the rovibrational level \v' — —15, J' — 1) 
as a function of the detuning, A„/j/, from this level (/ = 1 
W/cm^) at T = 2 ^K (solid red line) and T = 20 ^K (dashed 
blue line). 



For the strongly mixed levels, the photoassociation line 
shapes, shown in Fig. IHlfor \v' = —15, J' = 1), are thus 
governed by the natural width, about 4.5 MHz in Fig. [51 
and thermal broadening is of secondary importance even 
at a temperature of T = 20 fiK. Due to the relatively 
short lifetime of the level, the profile manifests also only 
a very weak asymmetry [44J. For regular levels such as 
v' = —6 the opposite holds, i.e., the thermal width is 
larger than the natural width. An increase in tempera- 
ture from 2 fiK to 20 fiK therefore has a noticable effect 
on the photoassociation rate, cf. Ref. [33| for a detailed 
analysis of the effect of thermal broadening on the peak 
rate coefficients. We conclude that photoassociation of 
strontium atoms into strongly perturbed levels, albeit 
challenging, is within reach for an experimental setup 
such as that of Ref. [l3| . 

After observing that photoassociation into resonantly 
perturbed levels such as v' = —15 or v' = —26 should 
be feasible experimentally, the transition moments from 
these levels into bound levels of the electronic ground 
state are examined in Fig. [T] Furthermore, Fig. |8] shows 
the modulus squared of the vibrationally averaged tran- 
sition moments governing the spontaneous emission co- 
efficients, cf. Eq. (fTO|) . and the branching ratios, cf. 
Eq. (fT4|) . While the level v' = —6 decays predomi- 
nantly, with a branching ratio of more than 80%, into 
v" = —3, a very weakly bound ground state level with a 
binding energy of 0.17 cm^^, the strongly perturbed lev- 
els w' = —15 or v' = —26 decay into a range of the ground 
state levels, including deeply bound ones. The largest 
transition moment is observed for the ground state level 
v" = 6 with a binding energy of 836.4 cm~^. The cor- 
responding branching ratios amount to about 17% for 
both v' = —15 and v' — —26, compared to less than 2% 
for v' — —6. Note that the branching ratios to v" ~ 6 
in Fig. m are almost equal for v' = —15 and v' = —26, 
while the transition moments in Fig. [7] are not. This is 
due to the dependence of the spontaneous emission co- 
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FIG. 7: (Color online) Modulus squared of the vibrationally 
averaged bound-to-bound electric transition dipole moments 
between excited state rovibrational levels v' — —6, v' — —15, 
v' = —26, all with J' — 1 (shown in Fig. [J) and all vibrational 
levels \v" , J" = 0) of the ground electronic state, X^E^. r 
denotes the lifetime for spontaneous decay to the X^E+ state. 
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FIG. 8: (Color online) Branching ratio for the spontaneous 
decay from levels v' = —6, u' = — 15 and v' = —26 to bound 
rovibrational levels of the ground electronic state. 



efficients on the transition frequency in addition to the 
transition moment, cf. Eq. ([TU|) . Based on the favorable 
transition moments between the strongly perturbed ex- 
cited state levels and v" — 6, stimulated emission using 
a nanosecond pulse could be employed in order to pump 
the excited state population selectively into the ground 
state level v" = 6. Alternatively, final state selectiv- 
ity could be achieved by photoassociation via Stimulated 
Raman Adiabatic Passage (STIRAP) [38]. It requires a 
sufficiently steep trap to ensure a well-defined phase of 
the initial state \E, J") which is expected to be feasible 
in a deep optical lattice |39|. Due to their large transi- 
tion moments for both pump and Stokes steps, the path- 
ways E ^ v' = — 15(— 26) -^ v" = 6 would be the most 
promising routes for STIRAP photoassociation from an 
optical lattice into deeply bound levels. 
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FIG. 9: Modulus squared of the vibrationally averaged 
bound-to-bound electric transition dipole moments between 
all rovibrational levels \v' , J' — 1) of the Ot potential and all 
vibrational levels \v" , J" = 0) of the ground electronic state, 
X^E^ (for other possible combinations of J' and J" the pat- 
tern is almost identical). 



A complete overview over transitions between the rovi- 
brational levels v' of the excited Oj states below the 
^S -I- -^Pi asymptote and all ground state levels v" is given 
by Fig. [3 For clarity, the figure has been separated into 
two parts, showing the highly excited state levels v' in the 
top panel and the lower excited state levels v' in the bot- 
tom panel of Fig. [51 Note that we find 110 excited state 
0+ levels v' below the ^S + •^Pi asymptote with J' = 1, 
i.e., v' = —6 corresponds to v' — 104, v' = —15 to 95, 
and v' = —26 to 84. Considering first levels close to 
the ^S 4- '^Pi dissociation limit, we notice that the last 
two excited state levels have extremely weak bound-to- 
bound transition moments. The next ten lower levels 
display a single peak in their transition moments, indi- 
cating pure Franck-Condon transitions close to the outer 
turning point. This is typical for weakly bound, regular 
levels. Transferring the molecular population to shorter 
bond lengths is extremely difficult for such levels and re- 
quires many excitation-deexcitation cycles [38| . 

The first strongly perturbed level, v' = —15 (or v' = 
95), leads to a prominent series of peaks in the squared 



transition moment matrix. Figure [9] indicates that also 
the neighbouring levels of f' — —15 are significantly 
perturbed. This would be important for pump-dump 
schemes using picosecond laser pulses [40|, |4l| . An ex- 
cited state wavepacket ideally suited for selective popu- 
lation transfer into v" — —6 is obtained by superimpos- 
ing levels v' = 92, . . . , 98. This translates into a spectral 
width of the photoassociation pulse of 15cm~^, corre- 
sponding to a transform-limited pulse duration of 1 ps. 
Note that a previous study considering only the exper- 
imentally observed weakly bound levels concluded that 
short-pulse pump-dump photoassociation near the inter- 
combination line transition is not viable [1^. The main 
obstacle is the quasi- i?~^ behavior of the excited state 
potential that leads to a reduced density of vibrational 
levels for very small photoassociation detunings. The 
number of vibrational levels present is then too small to 
obtain a truly non-stationary wavepacket [19j. However, 
the picture changes completely for more deeply bound 
excited state levels such as those around v' — 95. The 
spectral width of the pulse can easily be chosen such that 
several vibrational levels are within the photoassociation 
window, without exciting the atomic intercombination 
line transition that would lead to loss of atoms [41j . The 
advantage of a time-dependent photoassociation scheme 
in the presence of non-resonant coupling lies in the dy- 
namical interplay that arises between the interaction of 
the molecule with the laser light and the spin-orbit in- 
teraction. In such a situation, a dynamical enhancement 
of the final state population was found for strong dump 
pulses, indicating that the efficiency of population trans- 
fer is not determined by the transition matrix elements 
anymore |41| . 

A key question is how accurate our predictions are 
regarding the position of the perturbed levels such as 
v' = —15 or v' = —26. There is no doubt about the 
presence of such levels since it results from the cross- 
ing between the c'^n„ and A^Sj potential energy curves, 
and this crossing was confirmed by a recent experimen- 
tal study [2^. Our ab initio data reported in Ref. [2l| 
are able to reproduce the rovibrational energy levels 
for J' = 1 obtained from the fit of the experimental 
data to a Dunham type expansion [25| to within 0.64 
cm~^. Considering all experimentally observed levels 
with J' < 50, the root-mean-square deviation between 
theoretically calculated levels and the raw experimental 
data is 4.5 cm~^. Perhaps this value, ±4.5 cm~-^, should 
be considered as a very conservative estimate of the er- 
ror bars in the binding energies reported in the present 
study. The main sources of error in the binding energies 
are the inaccuracy of the c'^n„ potential and its spin- 
orbit correction, A{H). Scaling of the c^IIu potential or 
the A{R) coupling by ±5% leads to shifts in the bind- 
ing energies by 2cm^^ to 2.5 cm^^, in particular for the 
levels with strong singlet-triplet mixing. However, very 
good results for the A^S+ state and for the atomic spin- 
orbit splitting of the '^P and '^D multiplets, obtained in 
Ref. j21,], suggest the accuracy of the c^^H^ potential and 
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FIG. 10: (Color online) Vibrationally averaged bound-bound 
Raman transition moments as a function of the binding en- 
ergy of the intermediate Oj rovibrational levels for three dif- 
ferent pathways discussed in proposals for the measurement 
of the time variation of the electron-to-proton mass ratio, 
me/mp [10, U^- Note the different scales for the transition 
moments. 



the A{R) coupling to be better than 5%. Note that scal- 
ing the other spin-orbit couplings, ^i{R) and ^2{R), by 
±5% has a negligible effect on the position of the bound 
levels. This confirms our assessment of the estimated 
error bars of ±4.5 cm~^ as rather conservative. While 
such error bars might appear to be relatively large from 
an experimental perspective, they are not surprising for 
a system with 78 electrons, strong relativistic effects, and 
the A^I]+ potential as deep as 8433 cm^^ that arc found 
in the strontium dimer. 

The Franck-Condon parabola typical for transitions 
between regular vibrational levels [i^, |43] is absent in 
Fig. [5] This reflects the strong perturbation of the vi- 
brational spectrum of the excited state levels due to the 
spin-orbit interaction. A reasoning on possible optical 
pathways solely based on the shape of the adiabatic po- 
tentials will therefore give a wrong picture. To empha- 
size this point. Fig. [TOl presents the transition matrix el- 
ements for Raman transitions that are relevant in the 
proposal for the measurement of the time variation of 
the electron-to-proton mass ratio, m^/mp [lOJ, |ll|. The 
idea is to transfer molecules into the X^S+ ground vibra- 
tional level starting from the weakly bound ground state 
level v" = —3 (corresponding to v" = 60) that is pop- 
ulated by spontaneous decay from v' = —6, the lowest 
excited state level previously observed in a photoassoci- 
ation experiment [17[ • One could expect the efficiency of 



have large transition dipole moments with both v" — —3 
(or v" = 60) and v" = 0, yielding a high efficiency for 
Raman transfer directly from v" = —3 to v" = 0. The 
maximum Raman moments are found for u' — 14 and 
v' — 16, cf. Fig. |9l These levels are almost pure sin- 
glet rovibrational states belonging to the A^E+ potential, 
and are only marginally perturbed by the spin-orbit cou- 
pling. For all the pathways presented in Fig.[TU]the most 
favourable intermediate levels v' are those which are ener- 
getically the highest and yet almost unperturbed, i.e., the 
levels located just below the crossing between the c^n„ 
and A^S+ potential energy curves. This is easily ratio- 
nalized in terms of the strong transition dipole moment, 
dz ( A <— X) , of these levels and their relatively good over- 
lap with the rovibrational levels of the X^E+ potential. 
The decrease of the Raman transition moments for the 
deeply bound levels excited state levels, with v' < 10 and 
binding energies larger than 2000 cm~^, is due to shift of 
equilibrium positions of the A^S^ and X^S+ potential 
wells, cf. Fig. [H 

The Raman transition moments from w" = — 3 — > f ' = 



14/16 



— are larger than any of the moments 



for transfer from v" — —3 to v" — 27. Of course, even 
higher Raman transition moments are found for optical 
pathways to v" = that start in v" — 6, cf. bottom 
panel of Fig. (TUl the level that is populated by photoas- 
sociation into a strongly perturbed excited state level fol- 
lowed by spontaneous or stimulated emission or pump- 
dump photoassociation, as explained above. We thus 
conclude that a single Raman transition after photoas- 
sociation is sufficient to obtain molecules in the X^I]+ 
ground vibrational level. The least intensity of the Ra- 
man lasers is required for optical pathways starting from 
v" = 6, i.e., after photoassociation into strongly per- 
turbed levels such as v' = —15 or v' = —26. The path- 
ways starting from v" = 6 come with the additional 
advantage that the transition frequencies of the Raman 
lasers differ only by £"^"=0 — £t)"=6 ~ 225 cm^^ com- 
pared to 792cm-i for v" = 27 ^> i;" = or 1061 cm^^ 
for v" = —3 -^ v" = 0. We would like to stress here 
that all these conclusions concerning the Raman transi- 
tions should strictly be valid as the intermediate v' lev- 
els between 2000 cm~^ and 1400 cm~^ are located in the 
bottom of the A^S+ well where the potential is known 
precisely [2l|, [2^ , and are almost not perturbed by the 
spin-orbit interaction. This also means that doing high- 
precision Raman spectroscopy with these state should be 
feasible and the spectra will not be obscured by the spin- 
orbit perturbation effects. 



IV. SUMMARY AND CONCLUSIONS 



smaller than the efficiency of a two-Raman-step transfer 
employing an intermediate state, v" — 27. Inspection 
of Fig. [TU] reveals, however, that this expectation is not 
confirmed. There exist a few excited state levels, with 
binding energies between 2000 cm~^ and 1400 cm~^, that 



Based on state-of-the-art ah initio calculations, we 
have calculated photoassociation rates and spontaneous 
emission coefficients for the photoassociation of Sr2 
molecules near the ^S -I- '^Pi intercombination line tran- 
sition. We have also analysed bound-to-bound transition 
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moments as well as Raman transition moments connect- 
ing vibrational levels in the electronic ground state, rel- 
evant to achieve transfer into the X^S^ ground vibra- 
tional level. The vibrational spectrum of the coupled 
c^n„, A^E+, B^E+ excited state manifold is found to be 
strongly perturbed. Therefore, optical pathways cannot 
be predicted based on the turning points of the adiabatic 
potentials. Consequently, the theoretical analysis needs 
to fully account for the spin-orbit coupling of the elec- 
tronically excited states. 

For excited state binding energies of about 13 cm^^ 
and larger, up to 2000 cm~^, strongly perturbed vibra- 
tional levels are identified. The strong perturbations re- 
sult from the resonant interaction of the coupled vibra- 
tional ladders of the c'^n„ and A^I]+ states. For Sr2, 
these levels are found to be particularly well suited for 
the stabilization of photoassociated molecules to the elec- 
tronic ground state, either via spontaneous or stimulated 
emission. The photoassociation rate of the strongly per- 
turbed levels is calculated to be c omp arable to that of the 
lowest level previously observed [131 ^^ a temperature of 
T = 20 /iK and about one order of magnitude smaller at 
T ~ 2 /iK. We therefore conclude that photoassociation 
into strongly perturbed levels should be feasible with the 
currently available experimental techniques. 

Strongly perturbed levels display large bound-to- 
bound transition moments with deeply bound vibrational 
levels of the electronic ground state. If photoassocia- 
tion is followed by spontaneous emission, this will show 
up as a dominant decay into X^E+(w" = 6), although 
a large range of ground state vibrational levels will be 
populated as well. State selectivity of the ground state 
levels can be achieved by stimulated emission, either em- 
ploying STIRAP photoassociation in a deep optical lat- 
tice |39l| or pump-dump photoassociation with picosecond 
pulses |4dl4H. 

Identifying in the experiment the strongly perturbed 
levels of the c'^n^, A^I]+, B^I]+ manifold that are partic- 
ularly suitable for efficient stabilization to deeply bound 
ground state levels requires a spectroscopic search since 
even state-of-the-art ab inito methods cannot predict the 
positions of the rovibrational levels with precision better 
than a few wavenumbers for such a heavy system like Sr2. 



The theoretical precision is limited here mainly by uncer- 
tainty of the c'^n„ state and its relativistic correction, and 
can be reduced only after emergence of new experimental 
data concerning the c'^n„, A^S+, B-'^S+ manifold of Sr2. 

Finally, the crossing between A-'^S+ and c'^n„ poten- 
tials will be important not only for the initial formation 
of Sr2 molecules but also for any subsequent Raman-type 
transition proceeding via the coupled c'^IIu, A^I]+, B^Ej 
manifold of states. The presence of unperturbed lev- 
els of the A^E+ state, that are located just below the 
crossing with the c^n„ curve, leads to the somewhat 
unexpected result that the weakly bound X-'^S^ vibra- 
tional levels just below the dissociation limit show larger 
Raman transition moments with the ground vibrational 
level than with levels half-way down the ground state po- 
tential well. Direct Raman transitions to the ground vi- 
brational level thus become possible for both weakly and 
strongly bound levels. When utilizing these transitions 
for population transfer by STIRAP, deeply bound levels 
such as v" — 6 come with the advantage of a smaller 
frequency gap between the pump and Stokes pulse and 
significantly larger transition moments translating into 
lower pulse amplitudes. 

There are thus at least two good reasons for future 
experiments on the strontium dimer to employ strongly 
perturbed levels of the c'^fl^, A^E+, B^E+ manifold: ef- 
ficient stabilization to deeply bound ground state levels 
and large matrix elements for Raman transitions between 
ground state levels. Our calculations show these experi- 
ments to be feasible with currently available experimental 
technology. 



Acknowledgments 

We would like to thank Tanya Zelevinsky, Paul Juli- 
enne and Svetlana Kotochigova for many useful discus- 
sions. This study was supported by the Polish Ministry 
of Science and Higher Education through the project 
N N204 215539. Financial support from the Deutsche 
Forschungsgemeinschaft (Grant No. KO 2301/2) is also 
gratefully acknowledged. 



[1] A. D. Ludlow, T. Zelevinsky, G. K. Campbell, S. Blatt, 

M. M. Boyd, M. H. G. de Miranda, M. J. Martin, J. W. 

Thomsen, S. M. Foreman, J. Ye, et al., Science 319, 1805 

(2008). 
[2] M. Takamoto, F.-L. Hong, R. Higashi, and H. Katori, 

Nature 435, 321 (2005). 
[3] M. D. Swallows, M. Bishof, Y. Lin, S. Blatt, M. J. Martin, 

A. M. Key, and J. Ye, Science 331, 1043 (2011). 
[4] S. Falke, H. Schnatz, J. S. R. V. Winfred, T. Middel- 

mann, S. Vogt, S. Weyers, B. Lipphardt, G. Grosche, 

F. Riehle, U. Sterr, et al., Metrologia 48, 399 (2011). 
[5] A. Yamaguchi, N. Shiga, S. Nagano, Y. Li, H. Ishijima, 



H. Hachisu, M. Kumagai, and T. Ido, Appl. Phys. Ex- 
press 5, 022701 (2012). 
[6] P. G. Westergaard, J. Lodewyck, L. Lorini, A. Lecallier, 

E. A. Burt, M. Zawada, J. Millo, and P. Lemonde, Phys. 

Rev. Lett. 106, 210801 (2011). 
[7] M. Chalony, A. Kastberg, B. Klappauf, and 

D. Wilkowski, Phys. Rev. Lett. 107, 243002 (2011). 
[8] R. Ciurylo, E. Tiesinga, and P. S. Julienne, Phys. Rev. 

A 71, 030701 (2005). 
[9] S. Blatt, T. L. Nicholson, B. J. Bloom, J. R. Williams, 

J. W. Thomsen, P. S. Julienne, and J. Ye, Phys. Rev. 

Lett. 107, 073202 (2011). 



11 



[10] T. Zelevinsky, S. Kotochigova, and J. Ye, Phys. Rev. [27 

Lett. 100, 043201 (2008). [28 

[11] S. Kotochigova, T. Zelevinsky, and J. Ye, Phys. Rev. A 

79, 012504 (2009). [29 

[12] K. M. Jones, E. Tiesinga, P. D. Lett, and P. S. Juhenne, 

Rev. Mod. Phys. 78, 483 (2006). [30 

[13] A. Fioretti, D. Comparat, A. Crubelher, O. Dulieu, 

F. Masnou-Seeuws, and P. Fillet, Phys. Rev. Lett. 80, [31 

4402 (1998). 
[14] C. M. Dion, C. Drag, O. Dulieu, B. Laburthe Toha, [32 

F. Masnou-Seeuws, and P. Fillet, Phys. Rev. Lett. 86, 

2253 (2001). [33 

[15] S. B. Nagel, P. G. Mickelson, A. D. Saenz, Y. N. Mar- 
tinez, Y. C. Chen, T. C. Killian, P. Pellegrini, and [34 

R. Cote, Phys. Rev. Lett. 94, 083004 (2005). 
[16] P. G. Mickelson, Y. N. Martinez, A. D. Saenz, S. B. [35 

Nagel, Y. C. Chen, T. C. Killian, P. Pellegrini, and 

R. Cote, Phys. Rev. Lett. 95, 223002 (2005). 
[17] T. Zelevinsky, M. M. Boyd, A. D. Ludlow, T. Ido, J. Ye, [36 

R. Ciurylo, P. Naidon, and P. S. Julienne, Phys. Rev. 

Lett. 96, 203201 (2006). [37 

[18] Y. N. Martinez de Escobar, P. G. Mickelson, P. Pellegrini, 

S. B. Nagel, A. Traverso, M. Yan, R. Cote, and T. C. [38 

Killian, Phys. Rev. A 78, 062708 (2008). [39 

[19] C. P. Koch, Phys. Rev. A 78, 063411 (2008). 
[20] C. P. Koch and R. Moszyiiski, Phys. Rev. A 78, 043417 

(2008). [40 

[21] W. Skomorowski, F. Pawlowski, R. Moszynski, 

and C. P. Koch, submitted for publication (2012), [41 

arXiv:1203.xxxx. 
[22] H. A. Bethe and E. E. Salpeter, Quantum Mechanics of [42 

One- and Two-Electron Atoms (Academic Press, New 

York, 1957), p. 170. 
[23] B. Bussery-Honvault, J.-M. Launay, T. Korona, and [43 

R. Moszynski, J. Chem. Phys. 125, 114315 (2006). 
[24] A. Stein, H. Knockel, and E. Tiemann, Phys. Rev. A 78, 

042508 (2008). [44 

[25] A. Stein, H. Knockel, and E. Tiemann, Eur. Phys. J. D 

64, 227 (2011). 
[26] R. Napolitano, J. Weiner, C. J. Wilhams, and P. S. Juli- 
enne, Phys. Rev. Lett. 73, 1352 (1994). 



K. Sando and A. Dalgarno, Mol. Phys. 20, 103 (1971). 

B. Bussery-Honvault, J.-M. Launay, and R. Moszynski, 
Phys. Rev. A 72, 012702 (2005). 

T. Zelevinsky, S. Blatt, M. M. Boyd, G. K. Campbell, 
A. D. Ludlow, and J. Ye, ChemPhysChem 9, 375 (2008). 
V. Kokoouline, O. Dulieu, R. Kosloff, and F. Masnou- 
Seeuws, J. Chem. Phys. 110, 9865 (1999). 
K. Willner, O. Dulieu, and F. Masnou-Seeuws, J. Chem. 
Phys. 120, 548 (2004). 

S. Kallush and R. Kosloff, Chem. Phys. Lett. 433, 221 
(2006). 

C. Amiot, O. Dulieu, and J. Verges, Phys. Rev. Lett. 83, 
2316 (1999). 

S. Ghosal, R. J. Doyle, C. P. Koch, and J. M. Hutson, 
New J. Phys. 11, 055011 (2009). 

H. K. Pechkis, D. Wang, Y. Huang, E. E. Eyler, P. L. 
Gould, W. C. Stwalley, and C. P. Koch, Phys. Rev. A 
76, 022504 (2007). 

A. Fioretti, O. Dulieu, and C. Gabbanini, J. Phys. B 40, 
3283 (2007). 

R. Ciurylo, E. Tiesinga, S. Kotochigova, and P. S. Juli- 
enne, Phys. Rev. A 70, 062710 (2004). 
C. P. Koch and M. Shapiro, Chem. Rev. (2012). 
M. Tomza, F. Pawlowski, M. Jeziorska, C. P. Koch, 
and R. Moszynski, Phys. Chem. Chem. Phys. 13, 18893 
(2011). 

C. P. Koch, E. Luc-Koenig, and F. Masnou-Seeuws, 
Phys. Rev. A 73, 033408 (2006). 

C. P. Koch, R. Kosloff, and F. Masnou-Seeuws, Phys. 
Rev. A 73, 043409 (2006). 

M. Viteau, A. Chotia, M. AUegrini, N. Bouloufa, 
O. Dulieu, D. Comparat, and P. Fillet, Science 321, 232 
(2008). 

M. Viteau, A. Chotia, D. Sofikitis, M. AUegrini, 
N. Bouloufa, O. Dulieu, D. Comparat, and P. Fillet, 
Faraday Discuss. 142, 257 (2009). 

Note that the detuning in Fig. [6] is taken with respect to 
the binding energy _B„/^_i5 j/^^ = 12.9 cm~^, not with 
respect to the atomic transition frequency. 



